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Abstract. Thin films formed from small molecules rapidly gain importance in 
different technological fields. To explain their growth, methods developed for 
zero-dimensional atoms as the film forming particles axe applied. However, in 
organic thin film growth the dimensionality of the building blocks comes into 
play. Using the special case of the model molecule para-Sexiphenyl, we will 
emphasize the challenges that arise from the anisotropic and one— dimensional 
nature of building blocks. Differences or common features with other rodlike 
molecules will be discussed. The typical morphologies encountered for this 
group of molecules and the relevant growth modes will be investigated. Special 
attention is given to the transition between flat lying and upright orientation of the 
building blocks during nucleation. We will further discuss methods to control the 
molecular orientation and describe the involved diffusion processes qualitatively 
and quantitatively. 
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Figure 1. Examples of organic semiconductor molecules. In the first row the 
oligomers CgOi Porphyrin (H2P), and Alqs (tris(8— hydroxyquinoline)aluminium) 
are shown. In the bottom row examples for some rodlike oligomers are presented. 
From left to right: para— Sexiphenyl {6P), para— Quaterphenyl (4P), Anthracene 
(3A), Pentacene (5A), and Sexitiophene (6T). Please note that while 4P is shown 
in the twisted gas phase configuration, 6P is depicted fiat as it is found in the 
bulk crystal structure. 



1. Introduction 

Research over the last decades on growth processes at an atomic scale has greatly 
enhanced our understanding of thin film formation and crystal growth. In particular 
the realization that in addition to thermodynamic efi^ects also kinetic limitations at the 
surface play an important role during the growth of thin films and crystals has helped 
to explain many growth phenomena [IK5]. Although a great level of understanding has 
been reached for many different processes, the vast majority of the systems contained 
single atoms as the film forming entity. These are as such zero-dimensional particles. 

A new class of thin film materials — conjugated molecules (see figure [l] for 
examples) — has emerged in the past 20 years. As is demonstrated in figure [I] 
such molecules can be three-dimensional like tris(8-hydroxyquinoline)aluminium 
(Alq3), two-dimensional like the essentially planar Porphyrins and Phthalocyanins 
or one-dimensional like the Acences, oligo-Phenylenes and Tiophenes. The spherical 
Ceo molecule can either be viewed as a large zero-dimensional or a isotropic 
three-dimensional particle. Their use as organic semiconductors opens exiting 
new possibilities for electronic and optoelectronic devices, in particular flexible or 
stretchable [6] ones. There is ample of evidence that models used in inorganic epitaxy 
might also be applicable to this new class of materials fT'-'^ . However, on a molecular 
level these materials do not always follow the well established findings of classical 
epitaxy [101 E]. In this topical review, the underlying reasons and consequences for 
thin film growth will be discussed. A vast variety of different conjugated molecules are 
currently under investigation and their number keeps steadily increasing. As many of 
the effects described here are thought to be generic, we will limit this topical review 
to rodlike molecules. They are a representative group within the vast number of 
molecules with extended dimensionality. In fact, we will restrict ourself mostly to the 
oligophenylene molecule para-Sexiphenyl (6P) [T^HH]- This model molecule is widely 
investigated for its potential use in organic thin film transistors (OTFT), organic 
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Figure 2. Schematic comparison of required molecular orientations for device 
architectures based on rodlike conjugated molecules, using 6P as an example, (a) 
For OLED applications, the 7r-system should be oriented parallel to the substrate 
and the top electrode. For rod— or plate— like molecules, this usually requires a 
flat lying configuration, (b) In an OTFT, the conjugated molecules should be 
oriented in an upright orientation. This facilitates an isotropic electric transport 
from source to drain parallel to the gate electrode. For holes the requirements on 
molecular orientation are unchanged as only the direction but not the path of the 
charge transport changes. 

light emitting diodes (OLED) and solar cells. In these devices, a specific molecular 
orientation is essential for optimal performance. By comparison to other rodlike 
molecules, we will highlight the modifications that arise in the growth behaviour by 
changing molecular properties. 

After a brief introduction to small organic molecules, first the nucleation 
behaviour will be discussed. It is during this initial growth stage when the molecular 
orientation of the thin film is determined. This topic is of special importance as 
different applications require different molecular orientations on the substrate. The 
need to control this orientation is rooted in the anisotropic properties of the molecules. 
In particular, the charge carrier mobility for different crystallographic directions often 
shows a pronounced anisotropy |161 I17j . As a rule of thumb, charge transport is 
always best in directions with maximum tt/tt overlap. Often this is perpendicular to 
the long molecular axis or the molecular plane containing the conjugated 7r-system for 
two- and three-dimensional molecules. In figure[2j the required molecular orientations 
for a OTFT and an OLED are depicted for the model molecule 6P. As can be seen, 
the desired orientation for an OLED application would be fiat lying molecules with 
their 7r-systems parallel to the electrodes. Charge transport in such a device is then 
perpendicular to the 7r-system and the long molecular axis or largest plane. Also the 
desired light emission is maximised in such a configuration. On the other hand, for an 
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Figure 3. Schematic growth morphology for (a) VoUmer— Weber, (b) Frank— van 
der Merwe, and (c) Stranski-Krastanov growth modes. Reprinted with permission 
from [18]. 

OTFT application an upright standing molecular configuration is wanted to facilitate 
charge transport parallel to the substrate in an isotropic way. Therefore, possibilities 
to influence the nucleation behaviour and consequently the resulting film morphology 
as well as the molecular orientation will be discussed. 

Secondly, we will describe different growth modes and mechanisms observed 
during the formation of thicker films. Examples of Stranski-Krastanov and VoUmer- 
Weber growth modes as well as Frank-van der Merwe or Layer-by-Layer (LbL) growth 
will be presented [Udl]. Although difficult to achieve, the latter is often the desired 
growth mode to fabricate continuous films with homogeneous thickness. The so 
obtained smooth interfaces have a lower number of defects and yield generally a higher 
charge carrier mobility [19,- 22 . Which of the above mentioned thermodynamical 
growth modes is realized, depends on the ratio of the different surface free energies. 
The following requirement has to be fulfilled for any thin film growth to happen. The 
sum of the surface free energy tJi between substrate and adsorbate, and a between 
the adsorbate and the vapor are smaller then the surface free energy of the substrate 
CFs [m ■ Depending on the evolution of the change of the surface free energy 

Ao- = (T + (T., - as (1) 

during deposition, the three growth modes presented in figure [3] can be distinguished. 
For Act < at all times it is feasible for the substrate to be covered by the adsorbate 
layer. This growth mode is usually referred to as Frank-van der Merwe, Layer-by- 
Layer (LbL) or two-dimensional growth. In case that Act > at all times clustering 
will occur. This mode is called VoUmer- Weber, island or three-dimensional growth. 
In case that Act < for the initial deposited layers but changes to Act > for 
subsequently deposited material, the film will start to cluster after a thin uniform 
layer has been deposited. This growth mode is called Stranski-Krastanov growth 
mode. The latter is frequently found in inorganic heteroepitaxy for systems with a 
significant but not too large mismatch |23j . Unfortunately in organic thin film growth, 
the last two growth modes are much more common than the desired LbL mode. 

The above sketched thermodynamic description of thin film growth is not 
sufficient if the necessary diffusion processes are kinetically hindered. A more atomistic 
approach that includes the actual pathways is then needed to describe the observed 
morphologies accurately |24) . In particular, interlayer diffusion and the associated 
barriers play a decisive role. Ehrlich-Schwoebel or step edge barriers [23 1^ . 
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activation barriers for intralayer diffusion, and tlie anisotropy of these properties 
influence the final morphology. Depending on the absolute and relative sizes of these 
barriers different morphologies will be found in the resulting thick films. As a result of 
the sizable step edge barrier often found in these molecular films, growth phenomena 
like mound formation and rapid roughening are frequently observed. As the strong 
van der Waals interaction — typical for conjugated molecules — often dominates all 
other intermolecular and molecule-substrate interactions, three-dimensional growth 
is characteristic for organic semiconductor thin films. 

For all growth phenomena mentioned it is important to realize that already for 
one-dimensional molecules at least two different scenarios have to be distinguished. 
A smooth film grown in LbL mode might be useful for OLED applications when 
formed from flat lying molecules. On the other hand, if the molecules have an upright 
orientation, already a small number of layers grown in LbL mode at the gate dielectric 
will yield a decent performance in an OTFT configuration. The reason that a small 
number of layers will suffice is related to the fact that all important charge transport 
processes are confined to the first two monolayers [57]. 

While many studies focus on the submonolayer regime and interpret the behaviour 
of individual molecules, we follow a mesoscopic approach. This is justified by the fact 
that the behaviour of larger ensembles of molecules allows to infer information about 
the molecular level processes [TJ [S] . In addition, the investigated mesoscopic size range 
correlates well with the final device dimensions. The relevance of different changes 
in the properties at this mesoscopic length scale can thus directly be related to the 
device performance. 

1.1. Organic semiconductors 

The conductivity of organic crystals has been studied already in the early 
20*^ century [53] • However, only with the discovery of electroluminescene 
these materials received additional attention from the semiconductor research 
community [301 131] . Inspired by the Nobel prize awarded work of Alan Heeger, 
Hideki Shirakawa, and Alan MacDiarmid in the 1970s, many researchers focused 
on conjugated polymers which exhibit good conductivity if prepared properly [3^. 
In the 1980s, organic heteroj unctions and organic thin film transistors 34-3Sj 
have been demonstrated. The final breakthrough happened after the realization of 
high efficiency electroluminescence from organic light emitting diodes built both from 
polymers [37|[3S] and oligomers 1321 HO]- Nowadays, organic semiconductors are either 
already used or are about to enter the market soon in countless applications such as 
large-area-lightning, flexible solar cells, and displays. These devices are based on the 
integrated use of OLEDs, OTFTs, sensors, and organic photovoltaic cells [4TU45j . 

With respect to their growth behaviour, organic semiconductors are of interest 
for several reasons. As described above, classic surface science treats zero-dimensional 
particles with a few exceptions such as Si dimers [46'. The extended shape of the used 
molecules does not only allow them to obtain different orientations in space but also 
influences the way they interact with the surrounding. It is the extended electronic 
system that is responsible for the large intermolecular forces. The underlying van 
der Waals forces are small for the individual constituting atoms but can add up 
to a few eV for the entire molecule. The extended electronic system also helps to 
smoothen the effective corrugation of the substrate. The molecule averages over 
many possible atomic adsorption sites to find the molecular adsorption site with the 
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minimum energy. Directly related to this is a large number of internal vibrational 
degrees of freedom gTl 35]. These have to be considered when discussing the 
interaction at interfaces. These interfaces can occur between the condensed phase 
and the 2D or 3D gas phase, but can also be boundaries between different crystalline 
domains. Unfortunately, the large size of the building blocks results in a large number 
of translational domains. Together with the usually low symmetry of the crystal 
structure this can lead to additional disorder in the film, which in turn hampers 
the final device performance. For a general overview on the properties of organic 
semiconductors, the reader is referred to the following books [151 [50] . 

1.2. Rod-like conjugated molecules 

The three most important groups of rod-like conjugated molecules are the para-n- 
phenyls and the groups of acens and thiophenes (see lower part of figure [T]). The 
first two will be of particular interest here, as they show opposed properties in 
some important characteristics. The n-thiophene molecules are chiral and exhibit 
interesting chiral phenomena in thin layers but not in the bulk [51] . 

In this article, we will focus primarily on para-Sexiphenyl (6P) [T^ [5^H51] . It 
is important to realize at this point that the acenes and para-n-phenyls have a very 
different stiffness of their backbone. The single bonds in the n-phenyls allow a certain 
fiexibility of the backbone as compared to the acenes which possess a much stiffer 
backbone formed by two bonds. In the case of the n-phenyls, the phenyl rings can 
twist with respect to each other (shown for 4P in figure [I]) . This twist is observed 
for single molecules either on the surface or in the gas phase. However, in a bulk 
crystal the n-phenyls obtain a flat conflguration (as shown for 6P in figure [T]) [HI |55] . 
Typically, rod-like molecules like the acenes or phenylenes form a so called herringbone 
bulk structure. This packing motif is characterized by an alternating left and right tilt 
of the molecular plane around the long molecular axis. As a result, the long hydrogen 
terminated side of one molecules faces the flat side of the neighboring molecules where 
the TT-system is located. This conflguration balances the quadrupol moment of the 
molecules most effectively. It should be noted that this packing motif is very different 
for what is found for two-dimensional plate-like molecules that favour a planar bulk 
stacking. 

2. Forming a nucleus 

In this section, the important quantity of interest is the so called critical nucleus size 
i* . It is deflned as the biggest number of particles forming a cluster that will become 
stable by adding one more particle. While there are various way's to extract this 
number, the use of rate theories [5] in combination with scanning probe techniques 
has been proven to be extremely successful. The central result of the underlying theory 
can be summarized by the relation 



where A'^ is the island number density, F denotes the flux of incoming particles in 
numbers of particles deposited per unit time and surface area. The scaling exponent 
X = i* /{i* + 2) holds the dependence on the critical nucleus size. 




(2) 



v = voe 



Eo/ikBT) 



(3) 
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is the particle jump rate on the surface. Here, Eo is the energy barrier the particle has 
to overcome, T the temperature and ks the Boltzmann factor. The pre-exponential 
factor vq = 2kbT/h is often referred to as the attempt frequency. For all practical 
purposes, this is on the order of 10^^ in inorganic systems. As we will see later, 
for organic systems ly'o can deviate substantially from this value |3S]. 

A second approach is based on the scaling hypothesis stating that the island 
statistics in the steady state regime will depend — besides v and F — only on the 
coverage 8 via the mean island size [SS] . Using this assumption, it is possible to derive 
the critical nucleus size although this is not very easy to apply in practice. This result 
has been extended by introducing an additional scaling for the capture numbers |57| . 
However, most often the empirical scaling function of Amar and Family |58] is used. 
Recently a similar approach has been introduced by Pimpinelli and Einstein based not 
on the island size distribution but on the capture zone size distribution |59H61j . The 
presented scaling methods have the benefit that they require less data, and often a 
single experiment or even a single image can be sufficient to extract the critical nucleus 
size. A review and more in depth analysis of the strong points and weaknesses of the 
individual methods can be found in [1] [2] . 

The above considerations are only valid for low coverages after the initial transient 
nucleation regime but before the coalescence occurs. This intermediate regime is called 
the steady state nucleation regime. 

An additional obstacle — for all three methods — results from the fact that the 
molecules are anisotropic. As a result they can obtain different in plane but also out 
of plane orientations. As we will see later, the single shot methods also carry the risk 
to overlook interesting kinetic behaviour with respect to deposition rate or growth 
temperature. 



2.1. Obtaining the critical island size 

In the case of complete condensation, i.e., when re-evaporation can be excluded, the 
island number density N can be written as [5] 

TVcx (^^^0^"/'=-^°. (4) 

Here, Tjj refers to the temperature of the sample during growth — often called 
deposition temperature. The energy parameter 

En — —, (5) 

(i* +2) ^ ' 

can be split further into the activation barrier for diffusion E^ and the binding energy 
of the critical nucleus Ei. Provided sufficient data is available, a plot of In TV vs InF 
allows the extraction of the critical nucleus from the slope 

i* 

i* +2 ^ ^ 

Plotting the island density in an Arrhenius fashion as InA'^ vs. l/Tp, one can extract 
again from the slope 

i*E£) + Ei 
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information on the involved energies. An assessment of the quahty of the recorded 
data and the relevance of the extracted data can be obtained by comparing the results 
from the two different methods [51 IHj 

VoF +ap\nF ^ j/qt + (8) 

with yon and yoT being the y-intercepts of the two above mentioned plots. 

The above described rate equation based approach has been used extensively in 
the past to obtain information on the critical nucleus size. Typically sizes between 
2 and 4 are found for different rodlike molecules. For 5 A on Si02, a value of 3-4 is 
reported similar to 2-3 reported for 6P on disordered mica(OOOl) substrates |62) . 

However, due to the large experimental data set necessary, in many studies also scaling 
laws are used to determine the critical nucleus size. Most of the results are obtained 
by applying island size scaling [58]. However, capture zone scaling [59H6T] seems to 
provide more reliable results |64H66| for some cases. In agreement with rate theory, 
scaling laws typically yield values between 2 and 3 for 6P [SH |H] and somewhat 
higher values between 3 and 6 for 5A [SI [HI [H] on Si02 or cyclohexane terminated 
Si{001} [67]. However, in particular for the growth of 5A care has to be taken with 
respect to the applicability of these single shot methods because of fractal growth 
morphologies [55] . 

As will be demonstrated in the next section, the stability of the possible nuclei 
does not necessarily depend in a homologous way on the nucleus size. In fact, 
due to effective shielding some configurations can be more stable than others (see 
figure |5| . 5A and 6P are very similar regarding their herringbone packing in the 
bulk. Consequently, similar molecular configurations in the critical nucleus will have 
comparable stability with respect to other configurations. 

As can be seen from ([s]), in Q depends also on the attempt frequency Vq- 
From the y-intercept j/qt in the above mentioned Arrhenius plot, one can extract 
this quantity. It is important to note that for molecules this value does not always 
correspond to the one typically found in inorganic diffusion of zero-dimensional atoms 
(:^o = 1 X 10^'^ s~^). In fact, here the value of i^q can be much higher. Values up 
5.6 X 10^^ s~^ [69! are reported from thermal desorption spectroscopy (TDS). However, 
common values obtained from TDS (experimental and theoretical) are in the order of 
1 X 10^^ 47, 48, 69 74^ . Recent rate equation analysis yields a value of 2 x 10^^ for the 
diffusion of 6P on mica [62] . The explanation of these high values is given by transition 
state theory |75| . In this theory, the pre-exponential factor depends on the partition 
functions of the particle in the diffusive and adsorbed state. In contrast to an atom, a 
molecule possesses many vibrational and rotational degrees of freedom. In particular, 
the latter ones contribute only to the partition function of the diffusive state. As a 
result, the pre-exponential factor deviates from the well known 1 x 10^'^ s~^. The 
difference between values obtained by rate theory [55] and the sometimes extreme 
values obtained from TDS [51] arises from the different target phases. While in a 
typical growth experiment — used for the rate equation approach — the molecule stays 
on the surface in a 2D gas phase, in a TDS experiment the molecules enter the 3D 
vapour phase. However, the 3D phase has an even higher number of degrees of freedom 
compared to the 2D gas phase, resulting in different partition functions. 

Independent of the above considerations, an interesting question concerning 
islands formed by upright molecules remains unsolved at the moment. How — and 
at what point during nucleation and growth — do the molecules obtain an upright 



Nucleation and growth of thin films of rod-like conjugated molecules 



9 




^mm mwm 

Figure 4. 6P cluster binding energy on a 6P(001) plane as function of cluster size. 
Graphs for clusters formed by upright standing and lying molecules are presented. 
Snapshots of possible cluster geometries are shown for cluster sizes of 3, 5, and 
7 molecules. While the top row shows the clusters formed from upright standing 
molecules, the corresponding clusters formed by lying molecules can be found at 
the bottom. Please note that in the clusters formed by flat lying molecules all 
molecules have their 7r-system exposed. In the seven— molecule cluster — and to 
a certain extend in the cluster formed by 5 standing molecules — some molecules 
have their 7r-system saturated. |62l Copyright (2011) by The American Physical 
Society. 

orientation? By considering the binding energy of an upright standing 6P molecule to 
the 6P(001) plane of 0.21,eV[|]and comparing it to the one for a flat lying molecule of 
1.27 eV [75] it is plausible to assume that a single molecule will always obtain a flat 
lying configuration [77] . Eventually flat lying molecules will meet and form initially 
unstable dimers and trimers that will decay or continue to grow and become stable 
when big enough. At some point — to form a film of upright standing molecules — 
the molecules have to change from a lying to upright standing configuration. Recent 
molecular dynamics (MD) simulations illustrate the problem [62j . The graph presented 
in figure |4] shows the evolution of the cluster binding energy with increasing cluster size. 
The cluster binding energy is calculated by comparing a situation with i molecules 
in a cluster to the same amount of i molecules adsorbed in a lying configuration on 
the 6P(001) surface. When the energy difference between the cluster and the separate 
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Figure 5. Illustration of some of the possible configurations of the critical nuclei 
for a rodlike molecule in an upright configuration. Assuming a herringbone like 
configuration similar to the bulk structure, different possible nuclei configurations 
are sketched for i*=0 to i*=4. The yellow labels give the number of fully exposed 
TT-systems (red bars). Configurations with one and 5 molecules yield a minimal 
number of fully exposed 7r-systems. Reprinted with permission from |78| . 

molecules is negative, the cluster is themiodynamically stable. One can see that for 
i < 4 only flat lying clusters are stable, while clusters formed from upright standing 
molecules are entirely unstable. For cluster size between 4 and 14 molecules, clusters 
of upright standing molecules are stable but lying clusters would still be favored. 
Only for clusters bigger than 14 molecules 6P continues to grow in the required 
(GDI) orientation. However, the exact numbers will also depend on the details of 
the interaction with the substrate. In the MD simulation presented in figure |4j the 
clusters rested on the 6P(001) surface rather than on an amorphized mica substrate — 
discussed above — which is hard to simulate. In any case, the simulation results will 
be correct for second layer nucleation as well as all following layers. 

More information can be extracted from the graphs shown in figure [4] on the 
growth of standing molecules. Already for clusters of only two molecules {i* — 1) 
the energy gain is quite large (on the order of 10 x kT) and these clusters would be 
very stable. Furthermore, a simple geometric argument based on counting the number 
of fully exposed 7r-systems arrives at the same critical island size of i* ~ 3 ... 5 for 
upright standing molecules. This is also the value often reported in literature [521 - 
IM] , In figure [5] several possible configurations of critical nuclei from i*=0 to i*=4 
are sketched. The fully exposed 7r-systems are also marked and their number is 
given. A large number of exposed 7r-systems is energetically unfavourable, and the 
system would in general try to minimize their number. Besides the rough nature of 
the model it turns out that only for clusters of four {i* — 3) or more molecules the 
number of fully exposed tt -systems becomes smaller than the number of molecules. 
From this simple model we can conclude that once the balance between molecules and 
exposed 7r-systems swings towards the molecules the nuclei become stable. This is in 
good agreement with the above presented MD results. The consequence of this can 
be seen twice in figure [4] First, the clusters formed by upright standing molecules 
become stable around a size which allows for molecules that do not belong to the 
island rim. Second, the final slope of the two graphs is different. With increasing 
cluster size the number of energetically more favourable molecules which are not part 
of the rim increases faster for clusters composed from upright standing molecules 
than for the fiat lying nuclei. At least for the investigated size range, clusters of fiat 
lying molecules contain only molecules which expose at least one 7r-system to the 
vacuum. Nucleation processes involving small metastable clusters (e.g. non-epitaxial 
dimers) that eventually convert into larger stable clusters are also observed in inorganic 
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semiconductor epitaxy [79H8T] . 

The important question is, how do these flat lying molecules obtain an upright 
standing orientation. Besides the reorientation of the whole cluster other scenarios 
are in principle possible too. As we will see below, defects on the substrate surface 
influence the orientation of the molecules. First principle calculations reveal that 
already a single defect in an otherwise perfect surface can alter the molecular 
orientation from flat lying to upright standing ^82J. The initial flat lying nucleus 
can act in a similar way. Molecules arriving later will undergo a kind of defect 
nucleation and obtain an upright standing configuration. This is also observed in 
experiments. When the coexistence of structures formed from flat lying and upright 
standing molecules is observed, the first ones are often responsible for nucleating the 
latter [83l[84]. 

Furthermore, already small clusters formed from upright standing nuclei should 
grow faster than their counterpart formed from flat lying molecules. While molecules 
with all rotational orientations can be incorporated into a cluster of standing 
molecules, they will eventually have to rotate in the case of a cluster of flat lying 
molecules. In any case, the calculation only yields the energetically most favourable 
configuration. It does not contain information on the probability that it can actually 
form. Thus, although energetically favourable the lying nuclei is actually more difficult 
to form due to the rotational hindrance. However, this implies that no substantial 
energy barrier exists for the reorientation of molecules from flat lying to upright 
standing. Although no information is available on such a barrier, it can not be too 
large, since the final morphology is dominated by upright standing molecules. Such a 
morphology would be unlikely in a scenario where initially both orientations compete 
and a high barrier would exist to obtain the upright orientation. In such a competing 
scenario — where both types of nuclei can form — a cluster of upright molecules could be 
kinetically stabilized just because it can grow in size much easier. Such an attachment 
limited aggregation (ALA) for 6P has been observed and a critical nucleus size of 
i* = 7 ± 2 is found [5S] . This is in reasonable agreement with the numbers obtained 
in the above simulations for the transition from unstable to stable for clusters formed 
by upright standing molecules. In such a scenario, no reconfiguration of the cluster 
from lying to upright would be necessary. 

2.2. Tuning molecular orientation and the role of defects 

Despite the persisting problems with respect to nucleation, several groups have 
succeeded in controlling the nucleation behaviour. This is an important step 
towards realization of functional devices, since different functionality requires different 
molecular orientations (see section [l] and figure [2]). 

A convenient, but technically not very practicable, method to control molecular 
orientation is via the substrate. The observed changes in orientation go hand in hand 
with a change of the surface free energy of the substrate. Although the surface free 
energy is therefore the most obvious ordering parameter for the change from flat lying 
to upright standing, we will use a different approach here. It is easy to see that 
the crystalline structure of the substrate surface plays a crucial role for the in-plane 
molecular orientation. However, the degree of order present in the substrate surface 
can also be decisive with respect to an upright or flat lying molecular orientation. In 
fact, the effect of changing the surface structure can often dominate over the behaviour 
expected from a substrate according to its other physical properties. A model system 




Figure 6. AFM analysis of 6P on mica: (a) 3D representation of a 850 nm long 
chain of 6P crystallites on mica. The film has been grown in 35 s at a substrate 
temperature of 360 K. (b) OMBE grown 6P fibers on mica(OOOl). The film has 
a nominal thickness of 4nm, and the sample temperature during growth was 
held at 360 K. The fibers consist of long segments, (c) 6P islands on carbon 
covered mica grown at 330 K (nominal film thickness Inm). (d) 6P islands on 
sputtered mica grown at 330 K (nominal film thickness Inm). The insets in (c,d) 
are cross sections revealing a terrace height of 2.6 nm corresponding to the length 
of the molecule, (a); Reprinted from 1881 . With kind permission from Springer 
Science and Business Media. (c),(d): Reprinted from [92]. Copyright (2007), with 
permission from Elsevier. 

demonstrating this is the deposition of para-Sexiphenyl onto muscovite mica surfaces. 
Several groups have shown that for a large number of rodlike molecules — such as 6P, 
5A, and 6T — nanofibers formed by flat lying molecules are the dominating morphology 
on clean mica surfaces [Ml [5BH^ . 

In figure [6] typical morphologies obtained by Atomic Force Microscopy (AFM) of 
6P films grown on mica are presented. Figure [6|^a) represents a particular interesting 
case of nanofiber formation on mica(OOOl). Here, the spontaneous rearrangement of 
small crystallites on top of a wetting layer into chains of crystallites is observed [88] . 
We will discuss this interesting growth later in more detail. Changing the growth 
method from the used High Vacuum (HV) Hot Wall Epitaxy (HWE) to Organic 
Molecular Beam Epitaxy (OMBE) under Ultra High Vacuum (UHV) conditions does 
not not infiuence the existence of a wetting layer. However, as can be seen from 
figure |6](b) this change in growth conditions results in the formation of larger more 
uniform chain segments with a less pronounced internal structure |92H94j . The above 
mentioned rearrangement process happens in-situ during HWE growth. 

The influence of the vacuum conditions can also be seen in photoluminescence 
(PL) spectra obtained from these anisotropic flbers (figure [?]). In contrast to steady 
state PL spectra obtained from HV HWE grown 6P fibers [5S], the broad band 
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Figure 7. (a) Steady state photoluminescence of UHV OMBE grown 6P fibers 
on mica (0001). The broad band emission around 480 nm found in HWE grown 
fibers is missing. All curves are normalized to the 0-0 transition at 399 nm. PL 
measurements of OMBE grown films courtesy A. Kadashchuk (Data reproduced 
with permission from \95^.). 

emission around 480 nm is suppressed for UHV OMBE grown fibers. The missing 
band at 480 nm can be connected to defect states resulting from structural defects. 
The absence of this band suggests a higher quality of the OMBE grown films. 

The important conclusion of the above observation is the following. HV conditions 
are usually sufficient to obtain reproducible and well performing organic thin films. 
Given that the molecules are reasonable stable against oxidation and UV light 
the obtained structures do not change over time in ambient conditions. However, 
the morphology (and the resulting properties) of organic thin films can easily be 
influenced by other small molecules. Results obtained under (U)HV conditions should 
be carefully reviewed with respect to their validity under ambient conditions. In 
particular wetting layers and other possibly metastable structures far away from the 
bulk structure are sensitive to adsorbates. 

Figure |6]^c) and (d) show the result of two different approaches to reorient the 
molecules from flat lying to upright standing. The molecules forming the crystallites 
and fibers in figure [6j^a,b) have their long molecular axis parallel to the substrate. 
To obtain the morphology presented in figure a surfactanlj§] — namely carbon — 
has been predeposited on the clean mica(OOOl) surface. TDS has revealed that the 
full saturation coverage of surfactant completely suppresses the formation of a wetting 
layer [52] • In the presence of the surfactant, 6P films of upright standing molecules are 
formed. Similar experiments have been performed with 6P [6W and 4P films 
on gold surfaces. In both cases, a reorientation of the molecules has been observed. 
It has to be pointed out that prior to the flat-lying/upright transition the epitaxy 
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OMBE 293 K 
OMBE 80 K 
HWE 260 K 




§ The term surfactant is used here for a substance that infiuences the growth. 
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between the molecular film and the substrate is weakened. For the growth of 4P on 
Au{lll} this is very well documented. On a clean Au{lll} surface, the (211) contact 
plane of 4P has all molecules with their long axis parallel to the surface. Half of the 
molecules have their 7r-system parallel to the surface. After the addition of only 15% 
of a monolayer of carbon, the 4P contact plane changes to (201). The long axis of 
the molecules is still parallel to the substrate surface. However, in this 4P plane all 
molecules are tilted into an edge-on configuration, where the 7r-system is not fully 
facing the substrate anymore. At the same time, the in-plane molecules loose their 
alignment with the substrate. After a further increase of the surfactant coverage to 
0.5 ML, the molecules obtain an upright standing orientation. The new contact plane 
is the 4P(001). The weakened interaction with the substrate results in the formation 
of bent 4P nano-fibers with a (201) contact plane on gold and eventually the formation 
of mounds composed of upright standing molecules |961 197] . 

Finally, figure also shows 6P islands on mica formed by upright standing 
molecules. However, this morphology has been achieved by breaking the surface 
symmetry of the substrate by ion bombardment of the substrate |92j . The nominal 
1 nm thick film presented in figure |6jd) consists of islands formed by upright standing 
molecules as is demonstrated by the cross section presented in the inset. The 
same reorientation can also be achieved by heating mica substrates to elevated 
temperatures [23]. More insight into the possible root cause for the reorientation on 
mica has recently been obtained using TDS. Putsche et al. demonstrated recently that 
the existence of a wetting layer and the final orientation of the molecules is determined 
by the amount of potassium present on the surface [99]. The mechanism for the 
reorientation is similar to what has been shown above for a carbon predeposition. 

Although similar experiments exist for other small conjugated molecules, the 
above series of experiments is unique as it clearly shows the importance for a well 
ordered periodic substrate which can guide the molecules. The gradual loss of this 
guidance initially leads to an inferior quality of the crystals formed by flat lying 
molecules. Finally, after significantly disturbing the surface order by surfactants or 
sputtering, the molecules choose to form films where their long molecular axis is 
perpendicular to the substrate surface. The loss of order in the substrate can be 
the result of a surfactant layer, sputtering, or the molecules themselves acting as 
homo surfactants. The latter is in particular true if the molecules are adsorbed on the 
substrate in a hit-and-stick mode. The initial molecules can not align in an ordered 
manner and create a disordered substrate for the next layer [1001 1101) . As expected 
from this discussion, using an amorphous substrate such as Si02 will lead to the 
formation of films formed by upright standing molecules [Mlll02lll03j . However, often 
the layer of upright standing molecules itself can again act as a well ordered substrate 
and can promote the successive growth of structures formed by lying needles. See for 
example the needles forming on the 6P mounds grown on sputtered mica, as presented 
in figure [9];b,c) [75]. 

An extreme example of surfactant mediated growth of nanofibers [104| is 
presented in figure [Sj The mica surface has been pretreated with either water or 
methanol to remove the potassium atoms from the cleavage surface [1051 1106j . It has 
been shown by AFM that the so created rings are several hundred nanometer high 
and have a typical diameter of 4 jrm (see figure [sj^a)). The luminescence micrograph 
presented in figure [sjjb) and the corresponding simulation (figure [8](c)) reveal two 
things. First, the molecules are lying fiat on the surface, and second their in-plane 
orientation rotates along the ring. As only molecules having their long axis parallel 
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Figure 8. (A) AFM image (10 pm'^ by lOiim^) showing a 6P nano ring on 
pretreated mica. The inset presents the cross section along the indicated line. 
(B) Luminescence micrograph and (C) corresponding simulated image. In (D), 
the intensity along the rim (symbols) is plotted together with the theoretical 
predicted one according to Malus' law. Adapted with permission from |104| . 
Copyright (2003) American Chemical Society. 



to the polarization of the incoming Ught yield maximum fluorescence, one observes a 



3. Diffusion and thin film growth 

In the previous section we have discussed the initial stages of film growth i.e. the 
nucleation. After islands with a certain orientation have formed they need to 
grow, coalesce, and evolve into a complete film. Unfortunately, the asymmetric 
building blocks used, often favor the formation of rough and anisotropic surface 
structures [TBI I107H109| . Typically, one aims at a smooth film with a low number 
of grain boundaries and a flat surface. Grain boundaries represent inhomogeneities in 
the thin film and will have a negative effect on the transport properties |110H113] . On 
top of that they also affect the optical properties of the thin film devices such as nano 
fiber wave guides [1041 1114H116] . Different growth techniques have been employed 
to achieve a smaller number of grain boundaries. Notably the use of thin organic 
layers [15] or Self Assembled Monolayers (SAM) [117[ and Super Sonic Molecular 
Beam Deposition (SuMBD) [118| allow to influence the grain size without changing 
the substrate. 

As has been discussed in the introduction, there can be severe deviations from 
the thermodynamic growth modes (figure [s]), if the necessary diffusion mechanisms 
are kinetically hindered. The analysis of the morphology and if necessary structural 
information allows to extract important physical quantities such as step edge barriers 




Nucleation and growth of thin films of rod-like conjugated molecules 



16 



as well as parameters describing the growth mechanism. The generalized rules that 
can be deducted from the in-situ and ex-situ observations of film growth are valid for 
many molecules. 

3.1. Mound formation 

As we have discussed earlier, thin films of conjugated organic molecules which have 
their long axis perpendicular to the substrate are desired for OTFT like applications. 
Although several approaches exist to achieve this desired orientation, the resulting 
films are often characterized by the formation of growth mounds due to rapid 
roughening ^ [TTl (TH 11081 11191 1120] . The resulting rough interface is undesired as it 
negatively affects the charge carrier mobility [TMS^] . Nevertheless, these rough films 
have been studied in the past as they allow interesting insights into molecular diffusion 
processes and thin film growth mechanisms. 

The process of rapid roughening is a result of an imbalance in the surface diffusion 
over step edges. If the downward particle fiux over the step edge is dominating, the 
desired Layer~by-Layer growth mode is facilitated. However, in the case of a high 
step edge barrier for downward diffusion, which is described by a significant additional 
barrier — the so called Ehrlich Schwoebel barrier (ESB) [25ll26] — mound formation will 
occur and rough morphologies are the result. For the case of a very high ESB and 
realistic growth conditions, only atoms landing on the terrace h—1 will be incorporated 
into terrace h. The coverage Q of terrace h at time t for a given amount of deposited 
material h = Ft with the flux F can than be expressed by a Poisson distribution 

- '^^^ 77" 

e.W = l-e-^5:-. (9) 

Because Qi — 1 — e~'' will always be smaller than 1, the first layer — and also all 
subsequent ones — never closes. As a result steep trenches are observed. Such Poisson 
shaped mounds have been described in inorganic systems already early indirectly |121j 
and using real space methods |122| . An in-depth theoretical analysis has been 
conducted later by Elkinani and Villain who used the ancient Greek Zeno Paradox 
to describe a peculiarity of the mound formation |123l 1124] , In thin film growth, the 
Zeno Paradox describes a situation where narrow trenches between mounds get so 
narrow, that the probability for an atom or a molecule to land in the trench becomes 
increasingly smaller as the trench width decreases. After the unlikely event of a 
molecule entering the trench, the now narrower trench will have an even smaller 
probability to be filled by further molecules, and thus stays open. For this to be 
observable, a high step edge barrier has to prevent molecules landing on higher lying 
terraces from descending onto lower terraces deeper in the trench. As a result, the 
mounds get higher and higher but would never coalesce thus the substrate would not 
be covered completely. 

A detailed analysis of the involved atomic or molecular diffusion processes shows 
that in fact the time scales for the different basic processes play an important role [125) . 
First, there is the traversal time 

Ttr ~ Ajv (10) 

a particle needs to visit all the sites on an island. Here, A corresponds to the size of 
the island (measured in lattice sites). The fiux of incoming particles F determines the 
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Figure 9. para-Scxiphcnyl mound formation on sputtered mica. The evolution 
of the film morphology with growing film thickness (4 nm, 10 nm, and 30 nm) can 
be seen. AFM image size; Spm; z-scale: 20 nm, 35 nm, and 50 nm. From |7b| . 
Adapted with permission from AAAS. 



time 



At=l/{FA) 



(11) 



between the arrival of tlie particles forming the film. Finally, the residence time |125l - 

[SH] 



bL 



(12) 



describes the time a diffusing particle spends on an island with a characteristic size 
L {L is the island's circumference). Here, a and b are geometry dependent constants. 
While, the first term is on the order of the residence time (12), the second term 
accounts for the increase in residence time due to the step edge barrier. The ratio 
a = v' jv of the hopping rates for on terrace jumps ([s]) and 



-Es/ikgT) 



(13) 



for step edge crossings, can be used to obtain the additional step edge barrier AEes 
between the barrier for on island diffusion Eu and the barrier for step crossing Es- 
It is crucial to realize at this point that v' and therefore also A.Ees can only be 
effective values for interlayer mass transport. Different edge terminations can in fact 
have very different hopping rates and barriers. It is important to note that of all 
involved processes the one with the smallest energy barrier will be the dominating 
one, provided that the morphology connected with it is occurring frequently |24| . A 
very good review on mound formation can be found in ref. [2 . 

Prototypical examples of mound formation in thin films formed by rodlike 
molecules can be found in particular for pentacene [HI 1129] and para-Sexiphenyl [75] . 
Figure |9] gives an overview on the morphology of growth mounds formed by 6P on 
a sputtered mica surface at 300 K. Careful sputtering of the crystalline mica(OOOl) 
surface destroys the symmetry of the surface and results in an disordered surface 
with an unchanged chemical composition |76[ 192] . This modification reorients the 
otherwise flat lying molecules into an upright orientation. With increasing film 
thickness pronounced mounds start to form on the mica surface. As can be seen from 
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Figure 10. Individual hexagonal shaped 6P mound and corresponding cross 
section. Single 2.7 nm high terraces can be identified. 



figure [To] and the corresponding cross section the mounds are formed from upright 
standing molecules and exhibit an irregular hexagonal shape. The steep trenches and 
the change in curvature of the mound slope that characterizes this morphology are an 
experimental verification of the Zeno Paradox described by Elkinani and Villain [124) . 

An analysis of this growth behaviour |125l 11301 1131j allows to extract the step- 
edge or Ehrlich Schwoebel barrier active in such a system |76| . The step edge hopping 
rate in the system can be obtained from the top terrace diameter [125] 



I cx 



(14) 



The size of the top terrace is limited by the fact that for terraces larger than I 
nucleation will occur on top of it making it the second to top terrace. For the above 
presented film, I has been measured to be 40nm±20nm [7B]. The probability for such 
a nucleation event is related to the ratio between At ([Tl]) and the residence time r ( 12 ). 



The hopping rate for the on-terrace diffusion can be obtained from kinetic nucleation 
theory [5] which relates the nucleation density N — 1/A^ to the hopping rate ly and 
the flux F via The average island distance A is measured to be Ijim [75] for the 
film presented in figure ]9] The relevant time scales here are again At (11) and the 
traversal time Ttr ( [To] ). The fact that the so obtained value for AEes — 0-67 eV is 
30 times higher than the barrier for 6P diffusion on top of a 6P(001) terrace (0.02 eV 
determined by molecular dynamics calculations |76j ) gives rise to the pronounced 
mound formation in this system. Such high barriers are not uncommon for organic 
systems. Fendrich et al. report 0.78 eV for the Ehrlich-Schwoebel barrier of flat 
lying 3,4,9,10-perylene-tetracaboxylic-dianhydride (PTCDA) on PTCDA(102) using 
empirical potentials and the nudged elastic band method |132j . For both systems — 6P 
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Figure 11. Probability a for step edge crossing for different i* . Solutions for 
regime II (green line) are only valid for i* = 1. The curve exceeds the valid range 
(green hatched area) for a larger critical nucleus. Solutions for regime III (red 
hatched area) are well within the bounds r"**! <Si a <S. V^^l"^ . 



and PTCDA — the calculated diffusion barrier on the terraces is smaller by at least 
one order of magnitude. 

However, further careful analysis shows that the above used method is only valid 



for the case of i* = 1. As we have seen in section 2.1 (figure HI) this is rarely the case 



in organic systems [521 US] • Following the arguments of |126) we obtain 

Z-n'a'^' (15) 
with r = ly/F. Depending on the detailed balance of the above presented time scales 



( 10 1-( 12 ) one arrives at one of four possible regimes. For regime I where a <^ T ^ 



we obtain pure Poisson growth. Regime II extends from F ^ <C a ^ F , where 



Si — i*{2i* — l)/{2i*{i* + 1) + 2) < 1 for all i* . The exponents in (15) are than given 
by 

For larger values of a but still smaller than F^-^/^ the exponents take the form 

?■* 7* + 1 

For a > F"'^/^, we enter the regime of weak barriers and the mound cross section 
starts to deviate from the above presented wedding cake shape. Figure [TT] plots the 
evolution of a in regime II and III for different i* together with the extent of their 
validity. As one can immediately recognize, regime II is only valid for the case of 
i* = 1. The value of a w 2 x 10^^° can be expressed in terms of the step edge barrier 
using ([3]) and (13). The exponential prefactors vq and Vq are assumed to be equal. 



As the current regime II reproduces the scenario used in |76] we obtain a similar 
result of Ai?£5 — 0.58 eV. The difference is attributed to the factor of probability 
which has a small effect on the final result and is neglected for this overview. This 
result is only valid if the pair dissociation time Tdis ^ '^trl'^ ^ which is not the case for 
i* > \. We do find valid solutions for T^is ^ '''trh (this regime III in [126J and |133j 
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corresponds to regime I in |125j ). A smaller barrier of AEes = 0.36 eV is obtained 
in this regime. From figure [TT] one can see that there is a weak dependency of a 
on i* . However, the change is within the error bar of the experimental data used 
as input for the calculation. As discussed above, for values of a above or below the 
marked regions II and III one enters pure Poisson growth or the weak barrier regime. 
Please note that for larger i* , regime III splits into a fluctuation and a mean field 
regime with identical scaling exponents, and the evolution presented in figure [11] has 
to be carefully reviewed for i* > 2. The difference being that for a small number of 
involved particles the common mean field approach is not valid. The above presented 
method — taking into account the statistical nature of the initial nucleation — is needed 
in the case of small numbers. The complete phase diagram of — Ina/lnF vs i* can 
be found in |126j . However, taking into account the size and dimensionality of the 
building blocks we expect the fluctuation governed regime to be valid for larger i* 
values than for atomistic processes. 

However, the barrier — even for the same edge — must not necessarily be constant 
during thin film growth. Two things have been revealed for the initial growth of 6P 
on sputtered mica (see figure |6jd)). First, for the given film thickness too few second 
layer islands have nucleated. Second the island are only 2 nm high, indicating a larger 
tilt angle for the molecules. As a result a value of AEes = 0.26 eV is reported for 
the first layer ESB in [75] using the method presented in |125) . However following the 
arguments presented above, this calculation and the size of the resulting step edge 
barrier has to be carefully revisited. The analysis based on jl26| shows that only for 
regime HI a valid solution can be found. In regime HI the critical islands size for 
second layer nucleation 

T'^a" (18) 

can be calculated using the exponents 

7 = , ^ , and n=-—-. 19) 
I* + 5 I* +5 

An additional problem in this calculation arises from the fact that in 

L k + 2 

/ = 1 - ci^ (20) 

which relates the second layer island fraction / to the critical island size for nucleation 
we find the exponent k. This exponent k is known to be 5 for fiuctuation controlled 
nucleation in regime II in case that i* = 1 |125j . By comparing (15) with (18) in |126| 
and (21) and (22) in |125j we obtain the general expression for regime HI k = i* + 3. 



The result of these considerations used in (18), (19), and (20) is plotted in figure 12 



The red curve shows the expected a for different i* in regime HI. Although this curve 
ranges within its limits ax = a^/'^* ^ a ^ q,-i/3^-2/(3x) — q,^ fQj. small and 
large i* , the general condition for the validity of the fluctuation controlled regime 

X < (21) 

limits its extend to values of i* < 2. For larger values of i* the mean field approach 
from |133j becomes valid. Based on the expressions (8b) and (7) in |133j we obtain 



R. 



c2 



7-2 /I \ ' +1 



(22) 
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Figure 12. Probability for step edge crossing for different i* . The solution in 
regime III is plotted in red together with the range of applicability. For i* > 2 
the fluctuation controlled method becomes invalid and the mean field approach 
(blue) becomes valid. 



and 

(23) 

for the critical island radius (22) and the second layer island fraction (23). For the 
limit of high barriers — so that a <C 2/Rc — the exponent has the form m = 2i* + 6. 
The result is plotted in blue in figure [12] Based on this new analysis of the second 
layer nucleation we obtain a step edge barrier AEes ~ O.leV in the fluctuation 
controlled regime for i* < 2. For larger i* values, as we deal with here, we obtain 
AEes ~ 0.2 eV using the mean field approach. The latter result has only a very 
weak dependence on i* . The value of AEes for second layer nucleation of at least 
O.leV is substantially different from the above presented result obtained from the 
Zeno Paradox and the overall mound shape of thicker films {AEes — 0.36 eV). The 
observed change in molecular tilt angle for the first few layers [76] can explain this 
discrepancy. A molecule crossing the step edge with its long molecular axis roughly 
perpendicular, will bend over the edge while descending. As the molecular tilt in 
the (001) plane gets smaller the necessary bending and thus the required energy gets 
smaller, too. 

The situation for a nucleus formed by flat lying molecules is more complicated, 
since we can not estimate the order of the dimer dissociation time from the finite 
dimer energy. In addition sterical hindering during the nucleation starts to play a role 
and one enters the regime of ALA [85 . In principle, the above considerations can be 
extended into this regime. For such a situation the exponent in ^ obtains the form 
X = 2i* /{i* + 3) |134j . In addition, the growth laws — forming the foundation for the 
calculation — will play an important role. 

When the molecule crosses the step edge it performs a complicated sequence 
of twisting, rotating, and bending [77] ■ AH these processes are costly in terms of 
energy, and add to the final barrier height the molecule has to overcome. However, 
the molecule will take the pathway for which all contributions are added in such a 
way that the final barrier will be minimal. In particular, for a molecule crossing the 
step edge with the long axis roughly perpendicular to the edge, the bending energy 
is a significant contribution. A decreased tilt of the molecular backbone during the 
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crossing lowers this energy term and consequently the overall barrier height. 

The existing studies demonstrate that the prediction of barrier heights is 
complicated and full of pitfalls. For the presented good agreement has been 

achieved initially between theory and experiment. The analysis assumed a molecule 
that crosses the step edge with the long axis perpendicular to the edge [75]. However, 
after relaxing some of the constraints used, other trajectories — with lower barriers 
{{AEes = 0.34 eV)) — involving difficult Fosbury Flop-like movements at the step 
edge were found [77] . The above presented new analysis of the experimental data is 
in good agreement with the evolved molecular dynamics simulations that predict a 
complicated step edge crossing process and comparable barriers. Focused research in 
this direction is important, as mound formation and layer dependent ESB values are 
common in organic thin film growth. A level dependent ESB (e.g. as described above 
and in [751 1135j ) often goes hand in hand with a change in tilt angle of the molecular 
backbone. However, the proposed Fosbury Flop-like step edge crossing can not explain 
the experimentally observed step edge barrier reduction. Provided the ESB becomes 
small enough the initial layers can completely close |135j . Other examples include the 
growth of DIP on native SiO^, which is characterized by a transition from LbL growth 
to mound growth. This is explained by changes in the interlayer mass transport '119' . 

Care has to be taken when comparing experimental results with simulations and 
DFT based calculations. It is important to realize that the experimentally obtained 
barriers based on averaged mound shapes represent the effective barrier in the entire 
film. Using well defined step edges (the (100) in [TGllT^) allows detailed insight into 
the dynamics of the step edge crossing for the specific facet. However, the elongated 
hexagonal shape of the 6P islands (see figure 10 1 has at maximum 2 of these (100) 
step edges. In addition, these two (the (100) and the (100)) will also have different 
tilt directions (inward and outward tilt) with respect to the top surface, which in turn 
are both different from the tilt angle of the other unit cell facets (the vertical (010) 
facet) and other possible step edges. This and other peculiarities of molecular step 
edges (see also Appendix B in [77]) show how difficult it is to make precise predictions 
of experimentally obtained values for step edge barriers. 



3.2. Growth of three-dimensional islands and fibers 

Although smooth films are usually preferred, the crystalline and one-dimensional 
nano-fibers presented in figure [6]ja,b) are one out of many examples of a useful non- 
smooth morphology. It is important to remember that these anisotropic structures 
grow from flat lying molecules. The previous section dealt with upright standing 
molecules where no or only a weak anisotropy in the substrate plane can be expected. 
In particular, the fact that blue lasing 115. l.'?!) has been shown for these fibers and 
that they can be used as waveguides [137j opens several possibilities for applications. 

Two cases have to be separated here. While often three-dimensional fibers grow 
directly on the substrate (VoUmer- Weber growth), in particular the fibers found on 
mica grow on a metastable wetting layer (Stranski-Krastanov growth). In both cases 
the molecule-molecule interaction dominates over the molecule-substrate interaction. 
The difference between these two types of interactions is large enough to facilitate the 
rearrangement and reorientation of entire crystallites as entities. 

The rearrangement of crystallites containing more than 140000 molecules is 
observed during the HWE deposition of 6P onto crystalline mica(OOOl) at 360 K |88| . 
During the deposition of 6P, first a wetting layer is formed. With increasing coverage 
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Figure 13. AFM images of IIWE grown 6P fibers on mica(OOOf). (a) Fibers — 
formed by a rearrangement process of several individual crystallites — line up on a 
dislocation network (indicated by green dashed lines) present in the wetting layer, 
(b) Detail of 6P crystallite chains with a length of several jrm, grown on clean 
mica under HV conditions. The high-resolution inset reveals the presence of small 
crystallites decorating the fibers. Reprinted from |88| . With kind permission from 
Springer Science and Business Media. 



crystallites grow on this wetting layer. However, after a critical amount of 6P has 
been deposited a rearrangement takes place and fibers — formed by the already existing 
crystallites — become the dominating morphological feature. Figure [6]ja) shows such a 
chain of crystallites. The particular arrangement of the individual chains with respect 
to each other (see figure [TsJ^a)) is explained by a strain relaxation mechanism at work 
in this system f5Bl. The stress induced by the crystallites in the wetting layer leads 
to the formation of a defect network (indicated by green lines in figure [Tsja)) that 
guides the rearrangement process of the crystallites. It is important to realize at this 
point that during this rearrangement process the crystallites move as entities. This 
relocation of whole 6P crystallites on mica(OOOl) is possible due to the delicate balance 
between the strong intermolecular interaction and the rather weak film substrate 
interaction. Detailed x-ray diffraction (XRD) studies have revealed the epitaxial 
relationship between 6P and the mica(OOOl) surface [T^I138| . In particular, they have 
shown that once the formation of needles sets in, the initially compressed spacing of 
the (111) planes quickly relaxes towards the bulk value '139'. 

Recently, a bimodal size distribution for the crystallites on crystalline mica 
has been observed. However, this behaviour for ultra-thin layers at a slightly 
elevated temperature of 400 K is only observed after exposing the samples to ambient 
conditions. Using TDS before and after exposing the sample to ambient conditions as 
well as AFM revealed that the initial present wetting layer is transformed into small 
crystallites. This second generation of smaller crystallites forms between the already 
existing fibers or chains of crystallites. [140 For thicker films, the material from the 
wetting layer is most likely captured by the large number of existing big crystallites 
and fibers. In figure 13 b) a thick film where long 6P needles have formed is presented. 
Several small second generation crystals are visible in the inset of figure 13 [b). 

As we have seen for the case of crystalline mica versus sputtered mica, the 
substrate plays an important role in determination of the molecular orientation. 
However, also a particular surface reconstruction can provide an interesting growth 
template. The (1x1) reconstruction of the TiO2{110} surface is characterized 
by parallel rows of protruding oxygen atoms. These rows run along the [001] 
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Figure 14. 6P deposition on Ti02{110}-(1 X l).fT43l (a) High temperature 
deposition at 400 K leads to the formation of elongated islands formed by upright 
standing molecules, (b) Low temperature deposition at 300 K results in the 
formation of long 6P fibers. Please note that the orientation of the structures has 
rotated by 90° from [001] to [110]. (c,d) Cross-sections along the lines indicated 
in (a,b). 

azimuth |141| . The spacing of 6.5 A is sufficient to accommodate the width of a 
6P molecule. Deposition of 1.31VIL of 6P at 400 K leads to the formation of large 
islands formed from upright standing molecules presented in figure [l4| [ a). The islands 
are separated by trenches filled with small crystallites which are a few monolayers 
high. The trenches run parallel to the oxygen rows along the [001] azimuth of the 
TiO2{110} surface |142| . The islands themselves are polycrystalline with four domains 
symmetrically spaced around the [001] direction. With a size of only 300 nm by 30 nm, 
these domains are much smaller than the several jim large islands they form. The long 
axis of these domains is also oriented along the [001] direction |144j . The reason for 
this growth behaviour is rooted in the diffusion anisotropy present on this surface. 
The molecules can easily diffuse along the [001] direction guided by the oxygen rows. 
Analysing the width of the area in the trenches which is depleted from the small 
crystallites, one arrives at a ratio for the anisotropy between the diffusion along 
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(oxygen rows) 



Figure 15. Sticking anisotropy for 6P on Ti02- The fibers shown in the AFM 
image (z— scale: 30 nm) on the left grow along the [110]. The sketch shows the 
molecular arrangement from the top and the side. The grey arrow indicates the 
growth direction. Reprinted from |145| . Copyright (2006), with permission from 
Elsevier. 

[001] and [iTO] of 4 to 64 [US]. This is clearly a property of the TiO2{110}-(lxl) 
substrate surface since second and third layer islands show an isotropic shape. All these 
structures formed by upright standing molecules grow on top of a wetting layer of fiat 
lying molecules |146j . Such flat lying-to-upright transitions have been observed for 5A 
on Cu{110} [147] and other systems. For 5A on Cu{110}, the transition involves a fiat 
lying wetting layer, which is followed by an intermediate layer having a herringbone 
structure with the long molecular axis parallel to the substrate. For layers thicker 
than 2 nm an upright standing orientation is found in this system. 

Lowering the growth temperature to 300 K results in a complete change of 
growth morphology, molecular orientation, and mesoscopic structure orientation. The 
morphology presented in figure [T4{ ^b) is characterized by long and high polycrystalline 
6P fibers. It is important to realize that these fibers run parallel to the [110] and thus 
perpendicular to the oxygen rows and the trenches observed at higher temperatures. 
These fibers are formed from flat lying molecules that have their long axis roughly 
parallel to the substrate surface and are oriented along the [001] direction of the 
TiO2{110} surface. In addition to the diffusion anisotropy active at high temperatures 
here the sticking probability for molecules to be incorporated into existing fibers plays 
an important role. The long side walls of the fibers are terminated by the hydrogen 
atoms at the long end of the molecules. This has to be compared to the short side of 
the fiber where the tt systems of the molecules are exposed. It is clear that the sticking 
probability at the short end will be substantially higher. Consequently the fibers will 
grow quickly along [110] but slower in width. An illustration of the situation is shown 
in figure [TS} 

An interesting mesoscopic approach to orient the fiber growth has been shown 
by Madsen et al. |148j . They used arrays of gold coated micro-ridges. Tuning the 
ridge width and deposition temperature, 6P fibers growing perpendicular to the ridges 
could be grown with a high yield. 

3.3. Layer-by-layer growth of lying molecules 

Although the above mentioned needle-like morphology might be useful for special 
applications, a smooth interface is required for most applications. This is in particular 
related to the fact that a lower number of defects at the interface facilitates higher 
charge carrier mobilities |19fl22| . As we have seen above films formed by upright 
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Figure 16. Structural model of a 6P thin film grown on metal supported 
graphene (ligtit blue hexagonal layer) at 240 K. (a) The initial metastable layer 
is formed from flat lying molecules (grey carbon atoms) only, (b) After reaching 
a critical coverage, the structure changes to a bulk like molecular arrangement. 
This is achieved by tilting parts of the molecules on the long edge, as well as 
directly inserting molecules from the gas phase, (c) The film grows in LbL 
mode by a repetition of the previous two steps. The bulk of the thin film 
has a Baker-like structure 1151 and exposes the (111) plane to the underlying 
graphene substrate. The top most layer (orange carbon atoms) is not completed 
and shows a metastable structure consisting of flat lying molecules only. Adapted 
with permission from 151 . Copyright (2011) American Chemical Society. 



standing molecules suffer in many cases from high step edge barriers that will 
ultimately lead to mound formation and rough interfaces. So far, efforts to obtain 
Layer-by-Layer growth have led only to limited success. 

Wu et al. have achieved five layers in the desired LbL growth for the important 
case of 5 A on SiO^^ by using SuMBD |149j . Using conventional OMBD, Zhang et 
al. showed the strain relaxation driven transition from LbL to rapid roughening after 
5 layers for the plate-like molecule DIP |150| . In both studies, films formed from 
upright standing molecules have formed. However, some success has been obtained 
for flat lying molecules. 

Recently, Layer-by-Layer growth of flat lying molecules has been obtained for 6P 
on the technological important substrate graphene |151H153j . Graphene [1541 1155j can 
be used as a transparent, flexible and highly conductive electrode for organic electronic 
applications [1561 1157j . The combination of optical active films of flat lying molecules 
on a transparent electrode materials is a promising route to high efhciency OLEDs. 

The formation of the 6P film in LbL mode at 240 K has been monitored using 
in-situ real time Low Energy Electron Microscopy (LEEM) and micro Low Energy 
Electron Diffraction (yLEED) for structural characterization. The growth proceeds via 
a multi-step process that involves the reorientation of a significant portion of already 
deposited molecules |151[ . The process starts with the formation of a metastable 
layer of exclusively flat lying molecules (depicted in figure [16]^ a)). With ongoing 
deposition, this highly mobile initial layer [78l I153j transforms into a stable immobile 
layer having a higher packing density and a bulk like arrangement of the molecules 
(see figure 16 h)). The structure of this stable layer corresponds to the (111) plane 
of 6P. A similar growth process for the first monolayer has been reported for 6P on 
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Au{lll} t69j. However, on gold the growth at or above room temperature results 
in three-dimensional growth. It is important to understand the significance of the 
substrate for this process. Although earlier STM studies of 6P on graphite '158' also 
report flat lying molecules, the epitaxial relationship there is different to the one found 
on metal supported graphene. Using empirical force fields and total energy calculations 
it could be shown that indeed a different alignment of the long molecular axis is favored 
on the two substrates }151j . The Layer-by-Layer growth process continues with the 
repetition of the above two steps |151j . Every additional layer starts with the formation 
of the metastable initial layer of only fiat lying molecules that transforms into a layer 
with the bulk structure once a critical coverage is reached. Figure 16 c) shows the 
final film structure obtained by (pLEED) for 4.5 ML coverage. The achieved thin film 
structure bears the potential for high efficiency OLED structures on a transparent and 
fiexible substrate [151] . 



4. Conclusion 



In the first part, we have discussed the nucleation and growth behaviour of rodlike 
molecules. We presented several methods to determine the critical nucleus size. 
However, we also showed that due to the non-zero-dimensional nature of the 
molecules, care has to be taken when using formalisms originally introduced for atomic 
diffusion processes. Many problems — originally identified in inorganic systems — return 
in organic thin film growth. While phenomena like attachment limited aggregation 
have been observed in a few inorganic systems, they are encountered on a regular 
basis in organic systems. However, the biggest difference is due to the reorientation 
processes involved in the growth of films formed by upright standing molecules. Initial 
insight is gained mostly by using computational methods as the actual processes are 
difficult to monitor experimentally. The most difficult question to answer is related to 
the definition of the critical nucleus. Is an immobile cluster of flat lying molecules that 
finally nucleates a film formed by upright standing molecules the critical nucleus in a 
strict homoepitaxial sense? Although, this is to a certain level a semantic question, 
one has to realize that most experimental techniques to determine i* are insensitive 
to the orientation of the molecules. Consequently, this problem needs to be discussed 
when interpreting the results. 

Furthermore, the molecular orientation plays a crucial role in defining the 
efficiency of organic electronic devices. We have shown several ways to influence 
the orientation. We pointed out that substrate order and defects plays a crucial role 
for switching from films formed by flat lying molecules to films made from upright 
standing molecules. 

In the second part of this review we focused on diffusion processes that define 
the final film morphology. We extensively discussed mound formation in the presence 
of an effective Ehrlich-Schwoebel barrier and in particular the theoretical problems 
that arise when the critical nucleus becomes large. However, as expected a larger i* 
will not lead to a decrease of the involved step edge barriers. Two additional remarks 
have to be added to this discussion. First, given the size and the additional rotational 
degrees of freedom of molecules compared to atoms, the transition from the fluctuation 
determined regime to the mean fleld regime probably occurs at larger i* than in 
inorganic growth. Second, we briefly discussed how to extract the critical nucleus sizes 
for attachment limited aggregation (ALA) and diffusion limited aggregation (DLA). 
The latter is also observed in organic growth and leads to the formation of ramifled 
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islands. In the case of ALA, less molecules arrive at the lower step edge thus reducing 
the upward diffusion flux. For reasonable barrier heights this should facilitate LbL-like 
growth. In the case of DLA, the number of kink sites increases substantially, creating 
more low-barrier pathways over the step edge. While for ALA edge diffusion is in 
principle possible, this mechanism is not effective in DLA. Consequently, more and 
longer undesired domain boundaries are expected for DLA growth. The process of 
mound formation and its relation to Layer-by-Layer growth gets further complicated 
by the fact that the ESB can be layer dependent. In organic epitaxy often a change 
of molecular tilt angle is the root cause for this change of the barrier height and the 
growth mode. 

Self-organization of the formed nano-structures can effectively be controlled by 
balancing the anisotropics present within the growth system. This includes — but 
is not limited to — diffusion and sticking anisotropy as well as the anisotropy of the 
substrate on an atomic, but also mesoscopic length scale. Also wetting layers play an 
important role for many organic systems. They are not necessary stable under ambient 
conditions. However, rearrangement processes mediated by the wetting layer — like the 
one observed for the crystallite chain formation — can therefore only happen during 
(U)HV growth. No change of morphology by such a process is possible once the 
wetting layer has dissolved. 

In general, the anisotropy of the molecules leads to an anisotropy between the 
different diffusion processes. However, depending on the actual orientation of the 
molecules with respect to the substrate the diffusion pathways change their meaning. 
The path indicated by the red arrow in figure [iTjja) for the step edge crossing has 
a low probability and results in the undesired mound formation. However, for the 
case of flat lying molecules (figure [T7|^b)), the same diffusion process can actively be 
avoided. The small sticking probability at — what is now — the side of the fiber, allows 
the molecules to circumvent the red diffusion process by returning into the gas phase 
and reattachment at the small end. The consequence is the often observed fiber growth 
perpendicular to the long molecular axis. The same analogy holds for edge diffusion 
(upright standing molecules) and up hill diffusion (fiat lying molecules). 

Finally, we presented results of LbL growth of rodlike molecules and illustrated 
the often complicated rearrangement process occurring during growth of organic thin 
films (e.g. as for the above described low-temperature growth of 6P on graphene). 
Similar as above a deviation from the bulk structure goes hand in hand with a change 
of growth mode — in this case from three-dimensional needles to Layer-by-Layer growth 
of fiat lying molecules. The involved metastable islands exhibit an interesting diffusion 
behaviour which is mediated by a delicate interplay of strains in film and substrate. 

Considering the above presented information, two issues become immediately 
evident. Forcing the molecules to deviate from their desired bulk structure results in 
new and interesting growth phenomena and presents a viable route for controlling the 
film morphology. Secondly, the level of understanding of organic thin film epitaxy has 
increased dramatically within the last decade. However, there is still a large number of 
open questions. Having identified these questions, dedicated experiments — supported 
by computational methods — have to be designed to answer the existing challenges. 
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Figure 17. Sketch highlighting probable (green) and less probable (red) diffusion 
paths for rodlike organic molecules. The same diffusion paths with nearly identical 
barriers (ignoring the substrate effect) exist for structures formed from either 
upright or fiat lying molecules. While for upright standing molecules the red 
process is forbidden due to a high ESB, this process can be circumvented for 
structures formed by fiat lying molecules. 
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